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ABSTRACT 

We report the first direct detection of an exoplanet in the visible polarized light. The transiting planet 
HD 189733b is one of the very hot Jupiters with shortest periods and, thus, smallest orbits, which makes them 
ideal candidates for polarimetric detections. We obtained polarimetric measurements of HD 189733 in the B 
band well distributed over the orbital period and detected two polarization maxima near planetary elongations 
with a peak amplitude of ^2-10~ 4 . Assuming Rayleigh scattering, we estimated the effective size of the scat- 
tering atmosphere (Lambert sphere) to be 1.5±0.2^j, which is 30% larger than the radius of the opaque body 
previously inferred from transits. If the scattering matter fills the planetary Roche lobe, the lower limit of the 
geometrical albedo can be estimated as 0.14. The phase dependence of polarization indicates that the planetary 
orbit is oriented almost in a north-south direction with a longitude of ascending node £1=(16° or 196°)±8°. We 
obtain independent estimates of the orbit inclination /=98°±8°and eccentricity e=0.0 (with an uncertainty of 
0.05) which are in excellent agreement with values determined previously from transits and radial velocities. 
Our findings clearly demonstrate the power of polarimetry and open a new dimension in exploring exoplanetary 
atmospheres even for systems without transits. 

Subject headings: planetary systems — polarization — stars: individual (HD189733) 



1. INTRODUCTION 

More than 200 extrasolar planets have been discovered by 
the indirect methods of Doppler spectroscopy, photometric 
transi t s, microlensing , and pulsar timing (e.g. iMarcv et alj 
120051; lUdrv & Santos! 120061) . However, the direct detec- 
tion of exoplanets, enabling a study of their physical prop- 
erties, remains a challenge. Thus far atmospheres have 
been detected in transitin g planets (e.g. ICharbonneau et all 
l2002t iDeming et alJ 120061) and only in one non-transiting 
planet v And observe d in the far infrared with Spitzer 
dHarrington et alj|2006l). Polarimetry is praised as a power- 
fuLtechnique dHough et alj|2006t iKelleil (20061: ISchmid et alj 
2006) for detecting directly starlight that is scattered in a plan- 
etary atmosphere and, thus, possesses information on its ge- 
ometry, chemistry, and thermodynamics. However, due to the 
anticipated very low polarization degree of ^10~ 5 , a practical 
use of this method was generally attributed to the future. Here 
we report the first polarimetric detection of an exoplanet at a 
level exceeding theoretical expectations by an order of mag- 
nitude, primarily due to t he choice of the rec ently discovered 
hot Jupiter HD 189733b dBouchv et al.ll2005l) and the fact that 
the scattering planetary atmosphere is found to be somewhat 
larger than assumed before. 

The light scattered in the planetary atmosphere is linearly 
polarized perpendicular to the scattering plane. It is best char- 
acterized by the Stokes parameters q and m, normalized to the 
total flux. In general, when the planet revolves around the par- 
ent star, the scattering angle changes and the Stokes parame- 
ters vary. If the orbit is close to circular, two peaks per orbital 
period can be observed. The observed polarization variabil- 
ity should thus, in principle, exhibit the orbital period of the 
planet and reveal the inclination, eccentricity, and orientation 
of the orbit, and, if detected with high enough polarimetric ac- 
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curacy, also the nature of scattering particles in the planetary 
atmosphere. 

Thus, it appears feasible to detect the scattered light with 
polarimetric methods, even if an exoplanet is not spatially re- 
solved. Hot Jupiters orbiting their host stars on close orbits 
are therefore the best candidates for the first detection. Due to 
their proximity to the st ar, such planets apparently develop ex- 
tended hydrogen halos ( Vidal-Madjar et al. 2003), which may 
effectively scatter the stellar light, especially in blue wave- 
lengths. 

In this Letter we report the first direct detection of an ex- 
trasolar planet in the visible light. We present results of our 
polarimetric monitoring of the known exoplanet HD 189733b 
in the B band. A large number of measurements allows us 
to reduce the statistical error and for the first time detect 
directly the light scattered in the atmosphere of an extraso- 
lar planet. We interpret the data employing a model based 
on t he Lambert sphere approx imation and Rayleigh scatter- 
ing dFluri & Ber dyugina 2008) and deduce the radius of the 
sphere and the orbit orientation on the sky plane. 

2. OBSERVATIONS 

To increase the probability for detecting polarized scat- 
tered light from an exoplanet, we selected the nearby system 
HD 189733 with the planet at a short distance from the star 
(0.03 AU) and, thus, with a short orbital period (2.2 days). It 
was rec ently discovered in t he ELODIE search for transiting 
planets (Bouchv et al. 2005) and immediately became one of 
the favorite targets. 

Our observations were carried out i n 2006-2007 with the 
double image CCD polarimeter DIPol dPiirola et al.ll2.Q05h in- 
stalled on the remotely controlled 60 cm K VA telescope on La 
Palma, Spain. It is equipped with a rotating superachromatic 
half-wave plate as the retarder and a calcite plate as the ana- 
lyzer. To measure linear polarization, the retarder was rotated 
at 22.5° intervals. Each pair of the observed Stokes q and u 
was calculated from four exposures at different orientations of 
the retarder. To avoid systematic errors, the observations were 
carried out in cycles of 16 exposures, corresponding to a full 
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TABLE 1 

Parameters of the HD 189733 System. 
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FIG. 1 . — Polarimetric data (Stokes q and u with ±1<T error bars in the scale 
of 1CT 4 ) for HD 189733 and a Monte Carlo simulated sample. The original 
data are shown with filled circles: for HD 189733 in the top row panels for 
the year 2006 and in the middle row panels for 2007, and for the Monte Carlo 
simulated data combined for both years in single panels. The data rebinned 
for equal phase intervals are shown with open circles in separate panels. For 
HD 189733, the constant shifts in Stokes parameters Aq and An were sub- 
tracted from the data according to Table [T] The best-fit solutions deduced 
from the unbinned data are shown with solid curves. [See the electronic edi- 
tion of the Journal for a color version of this figure.} 

rotation of the retarder. In the 2006 season we made 10-15 s 
exposures at 2 x 16 positions, yielding eight single observa- 
tions of q and u per night. These were then averaged to calcu- 
late the nightly mean value and its standard error (la). Typ- 
ical errors of the 2006 measurements were 0.02-0.03%. In 
2007, in order to reduce the measurement errors, we increased 
the integration time for invidual exposures up to 20-30 s and 
made measurements at 4 x 16 positions, thus increasing the 
total integration time by a factor of 4. This reduced the errors 
by a factor of 2, down to 0.01-0.015%, which indicates that 
the accuracy was limited by the photon noise and did not suf- 
fer from systematic effects. Overall we obtained 93 nightly 
measurements for each Stokes parameter. This allowed us to 
reduce the statistical error down to 0.006% on average in the 
binned data (see Fig. [0 and clearly reveal polarization peaks 
of ~0.02% near elongations. 

For calibration of the polarization angle zero point we ob- 
served the highly polarized standard stars HD204827 and 
HD161056. To estimate the value of the instrumental po- 
larization, a nu mber of zero p olarized nearby (<25pc) stars 
from the list by iPiirolal (1 19771) were also observed. In fact, 
the instrumental polarization at the KVA telescope has been 
monitored since 2 004 within other projects as well (e.g., 
iPiirola et al.ll2005l) . These measurements demonstrated that 
in the B-passband it was well below 0.02% and invariable. 

3. MODELING 

To analyze the observed polarimetric signal, we employ a 
simple model based on the Lambert sphere approximation, 
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i.e., a perfectly reflecting surface with the geometrical albed o 
p =2/3, and Rayleigh scattering dFluri & Berdyu gina 2008). 
Modeling the observed variations in Stokes q and u allows us 
to reconstruct the orientation of the planetary orbit in space 
and estimate the effective size of the scattering atmosphere 
(Lambert sphere). In the model, fixed parameters are the or- 
bital period P, transit or periastron epoch 7q, semi-major axis 
a, and the radius of the star R*, which is considered to be a 
limb-darkened sphere. The values used are provided in Ta- 
ble [T| The limb-darkening was assumed according to IClaretl 
(2000), but its details were found to be insignificant within the 
measurement errors. Free parameters are the eccentricity e, 
orbit inclination i, longitude of the ascending node f2, radius 
of the Lambert sphere R^, and constant shifts in Stokes pa- 
rameters Aq and Am, which can be present in the data due to 
interstellar or circumstellar polarization. In the case of transit- 
ing planets, the orbit inclination can also be determined from 
photometric data, which is a valuable test for our model. Oth- 
erwise, polarimetry provides a unique opportunity to evaluate 
both i and fl Moreover, it is possible to distinguish between 
inclinations smaller and larger than 90°, which is not possible 
from transit data. 

In general, it is R^ that scales the amplitude of polariza- 
tion variations. The inclination scales the relative amplitudes 
in Stokes q and u. For example, at i = 0° q and u have 
the same amplitude. If /^O , the relative amplitude is also 
influenced by Q., e.g., at i = 90° variations appear only in 
Stokes q if fi = 0°, 90°, 180°, or 270°, and only in Stokes 
u if O = 45°, 135°, 225°, or 315° More examples can be 
found in iFluri & Berdvuginal ( 120081) . Observed polarization 
can be both positive and negative, since its direction is al- 
ways perpendicular to the line joining the planet and the star 
as projected on the sky plane. Our definition is in accordance 
with the common assumption that positive q is in the north- 
south direction, while the negative one in the east-west direc- 
tion. Positive and negative u are at an angle of 45° coun- 
terclockwise from the positive and negative q, respectively 
(see Fig. [3]). The inclination is defined in such a way that 
the planet revolves counterclockwise as projected on the sky 
for 0°< i <90° and clockwise for 90°< i <180°. Further, O 
varies from 0° to 360° starting from the north and increases 
via east, south, and west. 

In many cases, two maxima per period near the elongations 
are expected. For an orbit with moderate eccentricity the max- 
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FIG. 2.— The x 2 contours of the best-fit solution for HD189733b (solid 
lines). The contours are shown for the confidence levels 68.3%, 90.0%, and 
99.0%. Due to the intrinsic ambiguity in f! two equal minima separated by 
1 80° are seen. The Monte Carlo sample solutions are shown with dots. [See 
the electronic edition of the Journal for a color version of this figure.] 

ima shift closer to the periastron. Our polarimetric data reveal 
a clear periodic modulation with an amplitude of ^2-10~ 4 in 
Stokes q and about half of that in Stokes u, with the polariza- 
tion being negative, indicating the value to be near either 0° 
or 180° at high inclination, which is confirmed by the model. 

The model parameters are estimated using a \ 2 minimiza- 
tion procedure applied to the bulk of the original measure- 
ments for both years 2006 and 2007 simultaneously. Tests 
with simulated data showed that for fixed values of P and To 
and reasonably large number of measurements (>50) the pro- 
cedure is abl e to find a unique solution e ven for low signal-to- 
noise ratios dFluri & Berdvugi na 2008). The only ambiguity 
is in the f2 value of 180°, which is due to the intrinsic property 
of Stokes q and u to remain unchanged under the rotation by 
180°. The best-fit values of the parameters found at the nor- 
malized x 2 =104 are listed in Table Q] We find an excellent 
agreement with known values of eccentricity and inclination. 
Interestingly, the inclination of 98° revealed by polarimetry 
is very close to the complementary value 180°-86°=94° in- 
ferred with the transit method, which cannot distinguish be- 
tween inclinations smaller and larger than 90°. However, 
since the Stokes u amplitude is quite small, a better polari- 
metric accuracy is needed to distinguish between 86° and 94°. 
The fits to the observations and \ 2 contours of the solution are 
shown in Figs. Q] and [2] respectively. Note that the binned data 
also follow well the model, thus, indicating that the error dis- 
tribution is close to Gaussian and the binned data reveal the 
true variations. 

To evaluate the robustness of the solution to the measure- 
ment errors as well as to properly evaluate the errors of the 
model parameters, we employed the Mon te Carlo method. 
Following the recipe by Press et al] (119941) . we assumed the 
true orbit to be the same as our best-fit solution and simu- 
lated 200 samples of measurements (found to be enough for 
convergence) at the same orbital phases and with the same 
errors as the data for HD189733. One sample is shown in 
Fig.Q] We applied our \ 2 minimization procedure to the sim- 
ulated samples and obtained best-fit solutions for them. This 
allows us to obtain the distribution of deviations from the true 
set of parameters caused by measurement errors. In Fig. [2] 
the best-fit solutions for the simulated samples are plotted to- 
gether with the x 2 contours for the HD 189733 solution. Im- 
pressively enough the Monte Carlo tests concentrate near the 
X 2 minimum for all parameters, which proves that (1) the er- 
rors of our measurements have a Gaussian distribution, (2) 
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the signal is not spurious, and (3) the solution is robust to 
the errors. Thus, we are able to evaluate the uncertainties of 
the HD 189733 parameters using the Monte Carlo deviations. 
These are provided in Table Q] 

4. PROPERTIES OF HD189733B 

HD 189733b is known as a 'very hot Jupiter'. Its ra- 
dius in the B band and the orbit inclination w ere first deter- 
mined as fi=1.26±0.03flj and /=85.3°±0.1° dBouchv et all 
2005). Interestingly, the B band transit was observed 20% 
deeper than the transit in the R band, implying a larger ra- 
dius in blue wavel engths. Improved parameters were most re- 
cently obtained by lPont et all (120071) : #=1.154±0.017flj and 
/=85.68°±0.04°. However, the above estimates of the planet 
radius were obtained by measuring very accurately the planet- 
to-star diameter ratio and evaluating the stellar radius from 
the spectral class and colors with some u ncertainty. To over- 
come this uncertainty. iBaines et al.l(l2007l) measured the stel- 
lar disk directly using the CHARA Array in the near-infrared 
H band and yielded for the planet 1.19±0.08flj. The direct 
detection of the thermal emission with Spitzer using the sec- 
ondary eclip se indicated a 16jLtm brightness temperature of 
11 17±42K dDeming et alJl2006l) for the radius estimated by 
iBouchy et alJ ([2005). Very recently an infrared continuum 
spectrum with 0.49% ±0.02% of the flux of th e parent star 
dGrillmair et al.l2007l) and perhaps water vapour (Tine tti et all 
2007) were detected with Spitzer. 

The radius of the Lambert sphere /?l=1-5±0.2^j inferred 
by our modeling is ^30% larger than the radius of the opaque 
body revealed in the optical, although the lower ler margin is 
comparable with the latter and especially with the value by 
IBouchy etail d2005l) for the B band. Also, the Monte Carlo 
simulations indicate that the inferred value of might be bi- 
ased from the true value by ~0. 1-0.2 Rj due to an asymmetric 
shape of the minimum of the x 2 distribution (see Fig.|2|. The 
uncertainty in the radius is largely defined by the error in de- 
termining the polarization amplitude. In addition, since zero 
polarization at all orbital phases would correspond to the zero 
scattering radius (when the albedo value is not zero), mea- 
surement errors drag the solution to positive values and cause 
the asymmetry of the x 2 minimum. Nevertheless, the Lam- 
bert sphere radius corresponds to the largest possible geomet- 
rical albedo of 2/3, which implies that in reality the radius 
of the scattering atmosphere might be larger than 1.5,Rj. For 
instance, if the geometrical albedo of HD 189733b is compa- 
rable to those of the giant planets in the Solar system, e.g., 
0.52 for Jupiter and 0.47 for Saturn, then the scattering ra- 
dius of HD189733b would be R L =1.7±0.2Rj. A larger scat- 
tering atmosphere supports the idea of extended, evaporat- 
ing halos around hot Jupiters. If the ev aporation does take 
place , as in the case of, e.g., HD209458b dVidal-Madjar et alJ 

2003) , the size of the scattering surface might be compara- 
ble with the planetary Roche lobe size, Rrl- We estimate 
R RL for HD189733b to be 3.3 Rj dPaczvnski|[l97l and the 
corresponding lower limit of the geometrical albedo /?rl = 
2/3(/?l/^rl) =0.14. Another explanation for the detected 
polariz ation amplitude could b e a possible oblateness of the 
planet dSengupta & Maitill2006l) . 

An excellent agreement between the known and inferred 
values of the eccentricity and orbit inclination (Table Q]) indi- 
cates the plausibility of the assumed Rayleigh scattering, as 
its angular dependence strong ly influences the shape of the 
polar ization phase curve (e.g., Sea ger et"aT]|2000t IStam et ail 

2004) . Therefore, the scattering most probably occurs on 
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FIG. 3.— The orbit of HD189733b as projected on the sky. Solid and 
dashed lines indicate parts of the orbit in front of and behind the sky plane, 
respectively. The circle in the center depicts the star and the smaller circle 
on the orbit the planet. The direction of the orbital motion is indicated by the 
arrows. The reconstruction is made for i'=180°— 86°=94° and f2=16° (left) 
and 196° (right). If the inclination is 86°, the Q values are to be swaped 
and the direction of the orbital revolution is to be made counterclockwise. 
Positive directions of Stokes q and u and orientations on the sky are also 
shown. [See the electronic edition of the Journal for a color version of this 
figure.] 

small particles. Depending on the height in the planetary at- 
mosphere where the scattering is most efficient, it can be due 
to scattering on H, H2 or H2O, or even on small dust grains. 
If the latter is the case, the grain size should be <0.5fxm to 
have the Rayleigh type scattering in the blue. Such dust may 
be present in th e atmospheres of hot Ju piters, e.g., in form of 
silicate granes (Richardson et al. 2007). More accurate multi- 
color polarimetric measurements are needed to constrain the 
nature of scattering particles. 

Knowing all the orbital parameters of HD 189733b allows 
us to depict its orbit as projected on the sky plane and in- 
dicate the direction of the orbital motion. This is shown in 
Fig. [3] Such a plot is useful for future imaging or interfero- 
metric studies of the system as the coordinates of the planet 
can be predicted quite accurately. Since spectroscopy during 
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a transit revealed that the sky projections of the stellar spin 
axis and the orbi t normal are aligned to within a few degrees 
dWinn e t al. 2006), the direction of the stellar axis becomes 
known as well. 



5. CONCLUSIONS 

The first direct detection of the visible light from an extraso- 
lar planet allowed us to obtain fundamentally new information 
on the planet's orbit and scattering properties. We determined 
the orientation of the HD 189733b orbit projected on the sky, 
which can be used for subsequent direct detections by, e.g., 
radio interferometry. Moreover, our independent estimates of 
the orbit inclination and eccentricity are in excellent agree- 
ment with the values determined previously from transits and 
radial velocities. In addition, we inferred that the planet has 
an extended atmosphere which efficiently scatters the stellar 
light in the blue. 

Our findings open the door to new opportunities for direct 
detections of extrasolar planets, both hot Jupiters and Earth- 
like, to a large degree independently of their mass and gravita- 
tional effect on the host star. Furthermore, until now probing 
exoplanetary atmospheres was largely limited to systems with 
transits, which are relatively rare events. Polarimetry provides 
us with a new prospect to detect directly the light from the 
planetary atmosphere outside transits. Thus, with polarimet- 
ric accuracy approaching the photon noise limit, direct studies 
of exoplanetary atmospheres in the visible at any orbital incli- 
nation become reality. 
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